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Abstact: We investigate fluctuations in the momentum flux across a surface perpendicular to the 
velocity gradient in a stationary shear flow maintained by either thermostated deterministic or by 
stochastic boundary conditions. In the deterministic system the Gallavotti- Cohen (GC)relation for 
the probability of large deviations, which holds for the phase space volume contraction giving the 
Gibbs ensemble entropy production, never seems to hold for the flux which gives the hydrodynamic 
entropy production. In the stochastic case the GC relation is found to hold for the total flux, as 
predicted by extensions of the GG theorem but not for the flux across part of the surface. The latter 
appear to satisfy a modified GG relation. Similar results are obtained for the heat flux in a steady 
state produced by stochastic boundaries at different temperatures. 



1. Introduction 

There has been much effort during the past decade to connect statistical mechanics of stationary nonequilibrium 
states (SNS) with the theory of dissipative dynamical systems [1]. While most results obtained so far via this 
approach are more of a mathematical than physical interest there is one potential exception: the Gallavotti- 
Cohen theorem [2] and its generalizations [3] [4] [5] . The original GC theorem was motivated by the numerical 
results of [6] on a "thermostated" dynamical system. The phase-space time evolution of such a system is given 
by an equation of the form 

X^T{X) (1.1) 

with J- chosen to keep X{t) confined to a compact surface S in the phase-space while forcing the system into 
a nonequilibrium state. The latter requires that J- be non-hamiltonian with divJ-{X) = cr(X) ^ 0. 
Using some very strong assumptions on the dynamical system (1.1), GC proved that in the SRB measure [7] 
describing the SNS of this system the probability distribution Pt{p) — {S{p — TTr{X))) of 



1 

rja) 

satisfies the equality 



'^riX) = 373^ / cj{X{t))dt (1.2) 

/2 



]im ^\n^^ =p. (1.3) 

T^oo r((T) Fr(—p) 

Here (•) represents the average in the SNS. 

The quantity (a) is formally equal to the "rate of change" of the Gibbs entropy in the SNS. More precisely, 
if we start the system with a measure fj,o{dX) = po{X)dX where dX is the Liouville measure restricted to the 
surface S then, using the evolution (1.1), 

^G{t) = -H Pt\ogptdX = ^it{<j)^^{a) (1.4) 



dt _ 

where the existence of the limit will hold under the assumptions of the GC theorem. Furthermore we will have 
(a) < implying that Scit) —oo whenever the limiting state is not an equilibrium one with zero currents 
[8],[9]. 
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Based on the relation (1.4) ar{X) is often called the (normalized) "average entropy production during a time 
interval r" in the SNS. The identification of a^o, the object of the GC theorem, with entropy production was 
further strengthened by the form of a{X) in many examples of bulk thermostated systems, e.g. those considered 
by Moran and Hoover [10] for electrical conduction and by [11] for shear flow. In those systems <j{X) is given 
by an expression related to the hydrodynamic entropy production and the validity of the GC relation, eq.(1.3), 
was confirmed by numerical simulations, despite the fact that the conditions of the GC theorem are not satisfied 
there. 

Such bulk thermostated systems are however very different from realistic systems which are typically driven to 
SNS by inputs at their boundaries: the motion in their interiors being governed by Hamiltonian dynamics which 
do not produce any phase space volume contraction. It is therefore important, for comparison with real systems, 
to consider models of dynamical systems in which the thermostats forcing the system into SNS operate only 
near the boundaries. Such deterministic models were introduced by Chernov and Lebowitz [12] for shear flow 
and by Gallavotti [1], and van Beijeren [13] for heat flux. In this paper we investigate the validity and possible 
generalizations of the GC relation for such boundary driven systems. Before doing that we note that the GC 
theorem has been extended to open systems in contact with infinite thermal reservoirs which act on the system 
but are not changed by it. The simplest modeling of such a situation is via stochastic transitions, induced 
by the reservoirs, between different microstates of the system induced by the reservoirs [4] . Thus, to model a 
system carrying a heat current and/or a momentum flux, one may use Maxwellian boundary conditions. This 
means that a particle hitting the left (right) wall will be reflected with a Maxwellian distribution of velocities 
corresponding to temperatures {Tr) and mean velocities ul (ur) parallel to these walls. For Tl ^ Tr this 
will induce a SNS with a heat flux while ul ^ ur will (using periodic boundary conditions in the flow directions) 
induce a SNS with a shear flow, see [14] [15]. It is expected (hoped) that the deterministic (thermostated) and 
stochastic kinds of boundary modeling will yield similar SNS of a macroscopic system away from the boundaries. 
This is what happens in equilibrium systems, at least when not in a phase transition region. 

Such an "equivalence of ensembles" is far from established for SNS. In fact there is, at some level, a profound 
difference between thermostated and stochastically modeled SNS as far as Sq is concerned. As ready noted 
the former have Scit) -co, and Scit) (a) < 0, while the latter have Sait) — > Sg, Sait) 0. The 
origin of the difference lies in the differences in the measures describing these SNS. Thermostated SNS are 
described by an SRB measure which is singular with respect to the induced Lebesgue measure dX, while the 
SNS of stochastically driven systems are (this can be proven in some case and expected in general) described 
by measures that are absolutely continuous with respect to dX [16] [17]. This difference need not however mean 
much for a macroscopic system, since quantities of physical interest are sums of functions which depend only on 
a few variables. Their properties are therefore determined by the reduced distribution functions which can be 
expected to be absolutely continuous with respect to the local Lebesgue measure, i.e. expressible as densities, 
even when the full nic;asure is singular and fractal [18]. 

Interestingly enough, it is possible for the boundary driven pure heat flow case to model the thermostat in 
such a way that the expression for a{X) appearing in the GC relation is the same, up to terms whose average 
vanishes, for both the deterministic and stochastic case, see [1]. It can be written as 



where Jq{X) is the heat flux through some surface in the middle of the system and the average of the time 
derivative dF/dt vanishes in the SNS. (cr) has the form of entropy production in the left and right heat "reser- 
voirs" which is also equal to the hydrodynamic entropy production in the SNS [15]. 

The situation is different for the thermostated bomidary driven shear flow case considered by [12]. The (j(X) 
in (1.2), entering the GC theorem for that model is not clearly related to the hydrodynamic entropy production. 
The latter corresponds now to a momentum flux through the system, Jm{X), for which a GC relation holds 
for the stochastically driven system. The question is then whether there is still enough equivalence between 
systems driven deterministically and stochastically so that the GC relation for Jm{X), derived for the latter, 
also holds in the former. 

Another question which concerns both the heat conduction and the shear case is whether the GC relation can 
be observed in real macroscopic physical systems. More precisely, we know that in the linear regime the GC 




(1.5) 



so that 




(1.6) 
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relations imply an Onsagcr type reciprocity relations for the transport coefficients [19], [4] which hold not only 
globally but also locally. Our question is therefore the following: assuming that the GC relation holds for some 
flux crossing a surface S does it also hold (in some form) for the flux through part of S. The reason this is 
important for the applicability of the GC relation to real system is that for values of p for which Pr {p) is of order 
unity, Pr{—p): in eq. (1.3), which corresponds to the flux J going in the opposite direction from its usual one, 
e.g. the heat flowing from the cold to the hot reservoir, is so small in a macroscopic system that the possibility 
of observing it is effectively zero. A local flux reversal on the other hand may be quite observable: an attempt 
in this direction was indeed made by Giliberto and Laroche [20]. 

Here we describe numerical investigations of these questions for a deterministic and stochastically driven shear 
flow SNS and for a stochastic heat flow model. 



2. Description of the shear flow model 

Deterministic: The system consists of unit mass particles contained in a i x M box with periodic boundary 
conditions in the a;-direction and reflecting boundaries on the walls perpendicular to the y-direction. The 
dynamics in the bulk of the system is Hamiltonian with hard core interactions (the particles have a radius 
r). When a particle collides with the reflecting walls its outgoing speed is the same as the incoming one 
while the direction of the velocity is chosen in a way which simulates a moving boundary and creates a shear 
flow. The boundary transformation we consider is the same as in [12] [21]. Let ip and ■0 be the angle that 
the incoming (outgoing) velocity forms with the positive x-dircction if the particle collides with the upper wall 
or with the negative a;-direction if the particle collides with the lower wall. In the thermostated system the 
outgoing angle V is given by a function of the incoming angle (p: ip = f{ip). The function / we choose is 
f{ip) = (tt + 6) — -\/(7r + b)"^ — ip{ip — 2b). This is time-reversible, i.e. tt — /(tt — f{ip)) = if see Fig. 1. 



Fig. 1: Schematic representation of tiie dynamics of tiie system. 

Stochastic: The dynamics in the interior of the system is the same as before while the outgoing velocity at 
the upper (lower) wall is now chosen as if the particle was coming from a Maxwellian bath at temperature /3~^ 
moving at velocity vo{—vo). More precisely we assume that after a collision with the boundary the particle 
emerges with a velocity that is randomly chosen from a distribution: 



where Z is a normalization constant, Vq is the mean x-momentum of the particle after a collision and Cy is 1 if 
the collision is with the upper wall and —1 if it is with lower wall. 

In [21] we checked the validity of the GG relation for the phase space contraction generated during the collisions 
of particles with the deterministic thermostated boundary. We divided the phase space contraction into a 
contribution due to the lower boundary and one due to the upper one. We found that the GG relation was well 
verified for the total phase space contraction cr{X) but not for the partial ones^. 

We observe that based on the proof of GC we have no reason to expect that such a relation should hold for the partial phase space 
contraction. We tested it anyway since, as already noted, the GC relation appears to hold in more general situations than those covered 
by the GC theorem, e.g. for the total a here. 





(2.1) 
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As already noted there is no apparent connection between the phase space contraction and the hydrodynamic 
entropy production which is proportional to the flux of the x-component of the momentum across the system. In 
[12] equality between average phase space contraction rate and average hydrodynamical entropy production rate 
was shown to hold to first order in the shear in the limit in which the system becomes large (at constant density) , 
i.e. macroscopic, but the shear rate goes to zero in such a way as to maintain a constant total momentum transfer. 
This was done under the assumption that before a collision with either wall the particle velocities are distributed 
according to a Maxwellian. The equality between these averages was supported by numerical evidence. 

In the present paper we further investigate the possible equivalence between phase space contraction and 
entropy production. The momentum fltix is equal to the momenta carried by particles crossing a line through 
the middle of the system plus the exchange of momentum between two colliding particles when their center are 
on opposite side of the line, see fig. 2. 




Fig. 2: Events producing momentum flux. 



To be more precise we first introduce a discrete time for the system (in a slightly different way from what we 
did in [21]). Let X = {qi,Vi) a phase space point, ^t{X) be the time evolution induced by the dynamics and 
let t{X) be the first time at which a particle crosses the middle line or two particles on different sides of this 
line collide starting from the phase point X (we call such a situation a timing event) and let S{X) = $^(x)(-'^). 
Finally let tt{X) be the exchange of a;-momentum at the phase point X. 

We now specify the quantity whose fluctuations we will check. Given an integer r let 

r/2 

nAX)= '^(^'W) (2-2) 

i=-T/2 

and 



where the mean (•) is taken with respect to the stationary measure of the system. Let now n^(p) be the 
distribution function ofpriX) and 



We can formulate our "GC fluctuation relation for the entropy production" as follows: 



lim ^r{p) = Cp (2.4) 

T — ^OO 

Here C is the conversion constant between momentum flow and entropy production that, from hydrodynamics 

[21], is C = 2ui,/Th whc^re ut is the velocity of the particle; near the upper boundary and Ti, is the tempcratiire 
near the boundary. Observe that this as well as the following definitions, make sense also in the stochastic case 
when the phase space contraction rate is not defined. 

In this setting we define a local entropy production by looking at all events like the ones described in figure 2 
that happen in a specified part of the middle line of size IL. More precisely let 
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7r'(X)=7r(X)x[o,a]W (2-5) 

where xa{X) = 1 if the location of the momentum transfer event is at a point in A and otherwise. We now 
define in a way analogous to eqs.(2.2) and (2.3) the quantities ttI.{X), p\.{X), Il\-{p) and ^\.{p) and state our 
"local GC relation" for the entropy production as, 

lim^^(p)=Cp (2.6). 

r— >oo 

More generally we check if a relation of the form: 

\imerip) = Cip (2.7) 

r — ^cxD 

holds. 

The stochastic model permits us to discuss also different kind of transport phenomena. In fact wc can set the 
reciprocal temperature /? of the upper and lower walls to different value I3~^ and /3~ . In this case we will have 
also transport of heat through the middle of the system. We ran simulation also for this case, setting vq = 
for simplicity, and report the results in section 4. Clearly now the correct quantity to compute is the energy 
transfered across the middle line when a particle crosses or a collision happens. 



3. Numerical experiments for the shecir flow 

We simulated systems with N = 20, 40 and 60 particles and L, M such that the number density 6 = N/LM = 
0.034 was fixed. Wc chose; two different ways to set L and M . In one case wc fix L = Af i.e. a square domain. In 
the second case we keep M fixed at the value it had for = 20 and increase L proportionally as N is increased. 
In the deterministic case we fixed the energy per particle ^ J2i = 1 while in the stochastic case we fixed the 
values of vq and /? to reproduce the mean velocity Uh and temperature T5 = ((u — u)^) observed in the = 60 
simulation for the deterministic system. More precisely we fixed vq = 0.2 and (3~^ = 0.48. Finally the radius r 
was fixed to 1. For each value of N and L we followed a single trajectory of the system for 5 • 10^ timing events 
and used it to compute the distributions Il\.{p) ■ As in [21] (differently than in [22]) we did not discard any 
events between two consecutive segments of length t . 



3.1. Stochastic case 

As discussed in the introduction (see also section 5.1 for further discussion) we expect the fluctuation relation 
to hold for the total momentum flux corresponding to the hydrodynamic entropy production. This can indeed 
be seen in figure 3 in which ^r(p) is plotted for r = 300 and = 60 for the rectangular geometry. The dashed 
line represent the theoretical prediction ^oo{p) = 0.769p. Similar results hold for = 20 and A^ = 40. 

The relation in the strong form given by eq.(2.6) appears not to hold for Z < 1 but eq.(2.7) seems to hold as 
one observes a linear behavior of (,\-{p) always in Fig. 3. 

The results for can be used to obtain the behavior of the slope C; as a function of I. We report the result in 
Fig. 4 for the three valuc^ of A^ and the two different geometries we have used. Wc observe that in the case in 
which we keep the box square C'l decrease with A^ and seems to reach a limit different from 1 when N grows. 
If we just increase the horizontal side of the box keeping its height constant C; increases with N and it is not 
clear from the data what, if any, limit is reached when N ^ 00. 

Instead of fixing I we also tried fixing the length I * L but found nothing interesting. 



3.2. Deterministic case 

The situation looks very different in the deterministic case. No fluctuation relation seems to hold even when we 
consider the full momentum transfer. The result are shown in Fig. 5 where ^r{p) is plotted for r = 100, 200 and 
300. The dashed line represents the value predicted by eq. (2.4). Although we still observe a linear behavior, 
the slope appears to be increasing with r so that no limit seems to be reached. We will attempt an explanation 
of this phenomenon in the last section. 



5 



0.4 
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Fig. 3: Fluctuation relation for the momentum flux in the stoeastie shear flow with rectangular geometry. The filled circles (•) represent 
the experimental value for the total momentum flux, with error bars, for r — 400 and N — GO while the dashed line is the theoretical 
prediction. The other data represent the partial fluctuation relation for I = 0.9 (□), / = 0.6 (O) and / = 0.3 (X). In all three cases r = 400. 



We observe that a relation like eq. (2.7) seems to hold if we look at the partial momentum flux. This is clearly 
shown in Fig. 6 where the behavior of Cr(p) for I = 0.6 and several values of t is plotted for TV = 60 in the 
rectangular geometry. As in the stochastic case we can look at the behavior of C; as a function of I for both 
geometries. The results are plotted in Fig. 7. As can be seen the slope depends only very weakly on the size 
of the system, at least for the square geometry. The only value for which Ci seems to depend on the size is for 
I = 0.9. 

It is interesting to observe that if the fluctuation relation was true we would expect to observe a slope C; = 
2ub/Tb as discussed for the stochastic system. In this situation, and mainly for the square geometry, the value 
of Ub varies significatively from iV = 20 to = 60 while Ci remain almost constant. This and the result for 
the total momentum transfer suggests that the fluctuation of the phase space contraction rate and those of the 
momentum flux behave differently. 



4. Heat flow 

The stochastic boundary condition permits us to study the case in which the two walls are kept at different 
temperatures T+ and T_ . In this case the hydrodynamics entropy production is proportional to the heat current 
or energy flux from the upper wall of the system to the lower one. The events contributing to an exchange of 
energy are the same as those considered in Fig. 2 but now we consider the kinetic energy of a particle passing 
through the middle line or the exchange of energy in a collision between two particle that are on different side 
of the middle line. 

Analogously to what we did in section 2 we define s{X) as the energy exchange for a point X that is on the 
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Fig. 4: Slope Ci as a function of I in the stochastic shear flow for N = 20 (•), iV = 40 (X) and N = 60 (O). The left figure is for the 
square geometry while the right one is for the rectangular geometry. 
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Fig. 5; Fluctuation relation in the deterministic shear flow in the rectangular geometry for the total momentum flux with r = 100 (•), 
T = 200 (X) and t = 300 (O). 
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Fig. 6: Approach to a limit of the right hand side of eq.(2.7) for the deterministic case. In this case I = 0.6 while r = 100 (•), 300 (X), 
500 (O) and 700 (□). 




Fig. 7: Slope C; as a function of I in the deterministic shear flow for iV = 20 (•), iV = 40 (X) and N = 60 (O) in the rectangular 
geometry. The graph of the square geometry in analogous. 
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Poincare section with 



and 



r/2 
i=-r/2 



Let now -Br(e) be the distribution function of er{X) and 



Cr(e) = 7 — TVu 



As before we expect that 



lim C^(e) = Ce (4.3) 
where C is the proper conversion constant between heat flow and entropy production: 

c- — - — 

Similarly we define si{X), ei{X), -E|(e) and ^^(e) and check whether 

lim 4(e) = Qe (4.4) 

r — ^oo 

The numerical experiments are very similar to the ones described in section 3 but we considered only stochastic 
boundary conditions with rectangular geometry. Finally we fixed r+ = 0.7, T_ = 0.4 and all the other 
parameters as in the shear flow case. The global fluctuation relation is shown in figure 8. There it can be 
observed that the fluctuations are smaller than the shear flow case. In fact for r = 300, when in the shear flow 
case the limiting behavior was reached, we have just 2 points for the fluctuation. We interpret the results as 
showing an approach to the expected limit. Similar result are obtained for N = 40, 60. The plots for the ^(e), 
again, appear very linear, and their slopes C; are shown in figure 9. 



c, 
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Fig. 9: Slope Cj as a function of I in the heat flow model for AT = 20 (•), AT = 40 (X) and JV = 60 (<>). 
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Fig. 8: Fluctuation relation for the total energy flux in the heat flow model. The experimental value are plotted with error bar for 
T = 100 (•), T = 200 (X) and t = 300 (O) and N = 20. The dashed line is the theoretical prediction. 



Similar comments as the one for the rectangular geometry shear flow hold here. 



5. Conclusions 



We now try to summarize the results of our numerical experiments. 



5. 1 . Global fluctuation 



In the case of the stochastic boundary conditions the GC relation appears to be satisfied for both shear flow 
and heat conduction. This is not surprising. As in [4], [5], [15] we consider the stationary probability PT{X{t)) 
on the space of the trajectories of the system from time — t to time r. (This can be assumed to be unique, see 
for example [16], [23].) Define then the time reversal operation I{{qi,Vi)) = {qi, —Vi) then if X{t) is a possible 
trajectory so is I{X{t)) and we can consider the measure P = I*P. Then as in [4], [5] one finds 

^{X) = exp |i?(X(r)) - R{X{-r)) + a{X{t))dt^ (5.1) 

for some appropriate function (y{X). Here the Ihs represent a Radon-Nykodyn derivative and R{X) is a boundary 
term. The [2] relation then follows. 

It is easy to see that in the shear flow model Ut = J^^ ^^iX) is proportional to the momentum entering 
the system from the lower wall minus the momentum leaving the systc;ni from the upper wall. Due to the 
conservation of momentum in the bulk of the system the total momentum flux tTt through the middle of the 
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system is proportional to <Jr plus corrections due to the variation of momentum in the upper and lower half 
of the system. We expect these two quantities to fluctuate much less than the momentum flux, at least if the 
system is big enough so that we can expect to have a GC fluctuation relation for tt^ ^. Similar considerations 
hold for the heat conduction model. 

The deterministic case is far less clear. We know that the average phase-space volume contraction rate is equal 
to the hydrodynamic entropy production rate, only in some limiting situation, see [12]. Our numerical results, 
together with those in [21] show that this equality does not extend to their fluctuations. We observe, first of 
all, that although the CG fluctuation relation appears not to hold ^r(p) is linear in p so that we can rewrite 
eq.(2.4) as: 

^r{p) = CrP (5.2) 

The most striking effect appear to be the divergence of the slope Cr- A possible explanation of this can be 
based on the assumption that the distribution Yi-rip) to be close to a Gaussian, see [BCL], so that 

C 

where 

s{t) = ^j2 -^12 (5-3) 

t= — T t=—T 

is the integral of the 7r-autocorrelation function D(t) — (7r($'(-))7r(-)) — (tt)^. 

Now when X]t=-T ^i^) converges to a finite value the GC relation reduces to the usual Green-Kubo relation. 
On the other hand if we assume that ^(i) approaches when t ^ oo then Cr has to diverge as r^^. 

which is in good agreement with our numerical data. We were not able to check more directly this relation due 
to the lengthy simulations involved but we hope to come back to this in future work. 

We already noted that the quantity that enters the GC theorem for our deterministic shear flow model, i.e. the 
phase space contraction rate due to collisions, appears to have little in common with the hydrodynamic entropy 
production. This is different from models of deterministically boundary thermostated heat conduction systems 
in which the phase space contraction rate assume the form of an entropy production [1] [4], [13]. We expect that 
it is possible to introduce a deterministic forcing term for the shear flow acting only at, or near, the boundary 
and producing a phase space contraction rate equal to the hydrodynamic entropy production, see [12]. This 
would make the GC relation dependent on the form of the boundary forcing term. 



5.2. Local fluctuation 

As we already noted in [21] the local fluctuation, in both the stochastic and deterministic case, do not satisfy a 
fluctuation law in the form of eq.(2.6) but they appear to be in good agreement with the more general relation 
eq.(2.7). This seems to be contrary to the results obtained in [24] and in [5]. 

To resolve this apparent contradiction we note that in [24] Gallavotti considers a chain of weakly interacting 
Anosov dynamical systems^ and takes as a subsystem a finite piece of the chain. The phase space contraction 
rate is an extensive quantity as in the bulk thermostated systems. Furthermore the correlations in the chain 
decay exponentially both in space and in time. In a system with these characteristics one can prove that the 
fluctuation of the phase space contraction rate due to the degree of freedom of the subsystem satisfy a fluctuation 
relation with corrections proportional to the boimdary of the subsystem times the correlation length. 

In our situation none of the above characteristics is present. First of all we are not able to divide the degrees 
of freedom between the subsystem and the rest of the system. Moreover the phase space volume contraction is 
present only at the boundary of the system and the subsystem does not inchide any portion of the boundary. 
Finally we expect the correlation length in our system to be very long (potentially infinite), i.e. bigger that 
the size of the subsystem considered. For these reason we do not expect the arguments in [24] and [5] to be 



We observe however that differently from the analysis in [4] vjc do not have an a priori bound on the fluctuation of the total momentum 
of the system. 

The paper by Maes deals with a rather more general situation but similar arguments apply also there. 
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applicable to our situation. We observe however that our system is closer to the experimental situation described 
in [20] and to possible other experiments one can imagine doing. 

As in [21] we do not have any real explanation for the apparent validity of eq.(2.7). We think that the GC 
fluctuation relation can be extended to a partial relation of the form eq.(2.7) in a wide range of situations. 
These include a system of particles under the influence of an electric field and a Gaussian thermostat, see [22] 
and [25] for more details. We think that for the asymmetric simple exclusion process, one should be able to find 
an analytic justification of this behavior. 
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